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Abstract 

We introduce the concept of interpolation in quantum evolution and present 
a general framework to find the energy optimal Hamiltonian for a quantum 
system evolving among a given set of middle states using variational and geo- 
metric methods. The quantum brachistochrone problem is proved as a special 
case. 

1 Introduction 

Recently the Quantum Brachistochrone Problem (QBP) proposed by Carlini, et al. [1] 
has become a hot topicjl] The aim of QBP is to find the time optimal Hamiltonian 
under a given set of constraints for the quantum evolution between two given states. 

In this paper we will consider a more general problem: energy optimal interpola- 
tion in quantum evolution, and prove that QBP is a special case. 

The interpolation in quantum evolution can be described as follows: find a Hamil- 
tonian in a given Hilbert space under proper constraints, so that the quantum state 
lipit)) equals given states \ipi) ,■ ■ ■ , \4'm) at given times ti ,. . . , tm ■ The energy 
optimal interpolation (EOI) is to find the Hamiltonian in the solution space of the 
former problem such that 

tr{\H\^) = mm. (1) 

For the case m=2, the solution to EOI is the same as the result in If H is 
dependent on time, the evolution between ti and tj+i is reduced to the case m=2, then 
EOI turns trivial. In practical consideration, the quantum evolution is usually too 



^See, for example, [2,3,4,5,6]. 
^To be proved in section [31 
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quick to control the Hamiltonian for corresponding changes. Thus we only consider 
time-independent Hamiltonians. 

2 General Discussion 

Write the interpolation conditions in Schrodinger equation: 

exp{-iHti) lipi) = \ipi) ,i = 1, . . . ,m, (2) 

where we take ^ = 1 , ti = for convenience. 

A global phase factor will not alter the physical results, so ([2]) can be rewritten 

as: 

exp{-iHti) \iJi) = exp(i6'i) ,i = l,...,m, (3) 

where 9i can be arbitrarily adjusted to fit the physical system. 

Consider the Hilbert space of dimension n, in the eigenstate representation the 
Hamiltonian is diagonalizedl§ 

E = THT^ = diag{ei, . . . , (4) 

where T is the transformation matrix, T^T = TT"!" = I , ei ,. . . , En are the eigenvalues 
of energy. 

The state lipi) now transforms to \ipi) = e'^'T \ipi) . Write Anxm = dV'i') ; • • • ; li'm')) 
, then becomes: 

{Xkexp{-iekU))nxm = A (5) 
where (Ai, . . . , A„)"^ = T \ipi) , and normalization requires 

EiA.r = i (6) 

k=l 

Take the Hermitian adjoint of ([5]) and then multiply it by ([5]) we obtain: 

n 

(Y^ |A,pe-^(*-*^))^x™ = A^A = (e^(^^-^') (^#,))^x„^, (7) 
fc=i 

or writing in component form 

n 

^ |Afc|2 exp(iefc(ti - tj) + [{6^ - %)) = A^^-, l<j<i<m, (8) 

k=l 

■^In this paper we only concern Hermitian Hamiltonians. 
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where Aj^ = {ipi\ipj) = {ipj\ipi)* = A*^ ,so the case i < j is trivial. This is the 
fundamental equation for interpolation. 

From (IH]) we can immediately make the following observations. 

• The solution to the interpolation problem is invariant under a unitary transfor- 
mation applied to the given states. Hence two different sets of given states will 
yield the same result, if they are connected by a unitary transformation. 

• If (^i, . . . , Ai, . . . , A„, 6*1, ... , 9n) is a solution set to the equations, then for 
arbitrary Ae, {ei + Ae, . . . , £„ + Ae, Ai, . . . , A„, 9i — Aeti, . . . , 6'„ — Aet„) is also a 
solution set. The zero point of energy only cause a total phase change, thus we 
can always adjust the phases to adjust the zero point of energy without altering 
the solution. 

• There are m(m — 1) equations in (consider the real and imaginary part 
separately). The unknowns are 6i, |Afc| and e^. So there are in general — m 
equations and 2n + m unknowns in total. 

The structure of the solution space is determined by the number of independent 
equations. When the number of equations exceeds that of unknowns, we can use 
methods like least squares to find a path optimal Hamiltonian. When the number of 
unknowns exceeds that of equations, we can introduce the energy optimal Hamiltonian 
using constraint ([T]). 

As tr(|E|2) = tr{fH^f^fHf^) = tr{\H\y , (HD now becomes 

n 

J^4 = min. (9) 

k=l 

View the m(m — 1) equations in ([H]) as constraints, gives another constraint. 
Introduce corresponding Lagrangian multipliers aij,Pij,j, the variational function is 
defined as follows: 

n n 

k=l ^<j<i<''n k=l 
n n 

PijCZl l^'^l' sin(£fc(ti - t,) + 9i- 6,) - ^A,j)] + 7(5^ |Afcp - 1) (10) 

k=l k=l 

Variation with respect to Ek leads to: 



3 



dS 



l<j<i<m 



Pij cos(efc(ti - tj) + 9,- 9j)] = 0, /c = 1, . . . , n. (11) 

Variation with respect to |Aa;| leads to: 
dS 

— — =27|Afc| + 2|Afc| [aij cos{ek{ti - tj) + 9i - 9j)+ 

l<j<i<m 

Ai sin(£fc(ti - t,) + 9,- 9j)] = 0, k = l,...,n. (12) 

Variation with respect to 9i leads to: 
dS " 

^ = 5^ XI I^^H-^y sin(£fc(t* - tj) + - ^i) + 
* j^i k=i 

Pij cos(£fc(ti - tj) + 9i- 9j)] = 0, i = l,...,m. (13) 



where 



aji i< J, I « < J. 



Hence the fundamental equations for EOI are ([8]), (jlTl) . (1121) and (1131) . The un- 
knowns are Ek, \Xk\, (^i, ctij, Pij and 7. In total there are + 2n + 1 equations and 
rn? + 2n + 1 unknowns, and can be solved in principle. 



3 Quantum Brachistochrone Problem 

Constraint ([T]) is to minimize the module of H while the evolution time keep fixed. 
From Anandan-Aharonov relation which states that the 'speed' of quantum evo- 
lution is given by 2^H /h , Aif being the standard deviation of the Hamiltonian, we 
conclude the reversed constraint , viz. 

tr{\H\'^) = const. , t = min . , (15) 

should be equivalent to ([T]). 

Carlini, et al. ppj used the constraint 

tr{\H\'^) = tr{\H - tr{H)/n\'^) = const. ,t = min. , (16) 
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in solving QBP. The difference of H and H only lies in the choice of the zero point 
of energy. 

From (IH]) we know the the choice of the zero point of energy would not alter the 
solution, since tr{\H\'^) = tr{\H\'^) — \tr{H)\'^/n < tr(|H^p), constraint ([1]) requires 

tr{H) = 0, (17) 

i.e. tr(|ifp) = tr(|ffp). Hence ( IT5l) and (|T6l) are equivalent. So the the special case 
m = 2 for EOT should reduce to QBP. Now ([H]) and ffTTl) become: 

2ek =|Afcp t [a sm{ekt + 6) — [3 cos{ekt + 0)] , k = 1, . . . ,n; 

n (IS) 
J2\^k\'exp{i{ekt + 9)) = A, 

k=l 

for convenience we set ti = 6*1 = and omit redundant subscripts. Treat 6i^s and 
|Aa;|'s as known parameters we can solve (fT8l) . then adjust the values of ^j's and |Afc|'s 
to minimize the energy cost. 

To simplify f|T8|) . we shall now consider from the geometry viewpoint. The two 
given states and \ipf) span a space of dimension 2. It's natural that the time 
optimal evolution should be the geodesic in this space. So the Hamiltonian is of 
dimension 2 J| 

From (|T71) we have 

n 

J2ek = tr{E) = tr{H) = 0. (19) 

k=l 

Thus 

ei = -e2 = e. Sk = 0,k = 3, . . . ,n. (20) 

Then (ITSl) reduce to: 

2e = \Xi\'^t[asm{et + e) - l3 cos{et + 6)]; 
-2e = \X2\^t[asm{-£t + e) - (3 cos{-£t + 6)]; 

0= \\k\'^t[asme - pcose],k = 3,...,n; (21) 



n 

Ae-'^ = |Ai|V^* + IA2I V'^* + 5^ lAfci 



k=3 



'^See [8] for a proof. 

^We suppose the optimal Hamiltonian for m given independent states is of dimension m, which 
is a unproved guess. 
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The third eqn. of (pij) yields |Afc| = 0,k > 3. Substitute it into the forth eqn. we 
have: 

Ae'"* + (1 - A)e-'"* = Ae''^ (22) 

where we have used and write |Aip = A for convenience. 
Take module on both sides of ( 122|) : 

et|Ap = 1 - 4A(1 - A) sin(et)^ (23) 
Thus e = arcsin ^ 4A(1-a) "/'^' minimize e, we take A = 1/2. So 

e = arcsin a/1 — jAp/t = arccos |A|/t, (24) 

n 

tr(|^n = 5^ 4 = 2(arccos \A\/tf. (25) 

A:=l 

In the eigenstate representation, H = diag{e, —e, 0, . . . , 0), = (Ai, A2, 0, . . . , 0)"^, 
\'ip2) = exp{—iHt) lipl) = (Ai exp(— iet), A2exp(i£:t), 0, . . . , 0)"^. The Gram-Schmidt or- 
thonormalized state \ip2) = (|'^2) — coset sinet = (— iAi, iA2, 0, . . . , 0)"^, so 

(V^il - 1^1) (^2!) = diag{e, -e,0,...,0) = H. (26) 

At arbitrary time r, 

\^{t)) = (Ai exp(— ier), A2exp(i£:r), 0, . . . , 0)^ = coser + sin err \ip2) . (27) 

f l25p . fl26l) and 1^7} are the same as the main results in [T] which take tr{\H\'^)/2 = u"^. 
Substitute the results back into (1^ we have: 

|A| = Ae-'^ (28) 

thus {ip2\4'i) = Ae~'^ is real. The global phases of ipi and ip2 are adjusted to make 
(^2 1^1) be a real number at optimal solution. 

4 Case Aij = Sij, ti = {i — l)t 

Now we consider the special case when the given states are orthogonal and the given 
times are evenly spaced, viz. 

Ai, = 5ij ,U = {i- l)t. (29) 
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Since the given \ipi)^s are lineal independent. We have: 

m < n (30) 

Now ([H]) reduces to: 

n 

5^ lAfcP exp(i(fefet + 61)) = 0, l = t-j = l,2,...,m-l. (31) 

k=l 

In section [3] we concluded that the global phases 9i are adjusted to make Aij real 
at optimal solution. This conclusion is supposed to be still valid here. Thus we have 
01 = 0. 

When / = 1 a special solution can be easily found: the n vectors exp{iekt){k = 
1, . . . ,n) distribute uniformly on the unit circle, i.e. 

2k7c 

£kt= + ^0, lAfcP = 1/n, (32) 

n 

If n is prime, for any / < n, exp{ilekt) are still n vectors distributing uniformly on 
the unit circle. If n is not prime, when 1 is a factor of n, exp{ilekt) are n/l vectors 
distributing uniformly on the unit circle. In both cases we still have (1311) . But if 
/ = n, ( pTl) no longer holds, this is consistent with ( |30l) . 

The result is a special solution to the quantum interpolation, whether it is the 
solution to EOI depends on whether it is consistent with ffTTl) . ffT^ and ffTSj) . which 
remains unproved. 



5 Conclusion and Discussion 

We have educed the fundamental equations for the quantum interpolation and EOI. 
The general behavior of its solutions are preliminarily studied. QBP is discussed 
in detail as a special case. Another simple case for quantum interpolation is also 
considered. 

As future developments, the fundamental equations for EOI remains to be thor- 
oughly investigated, and this framework can be extended to mixed states and non- 
Hermitian Hamiltonians. The behavior of path optimal HamiltoniarJ^ is also worth 
studying. Meanwhile applications in quantum computation are probable. 

^See pagelH probably defined via methods such as least squares. 
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